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. Abstract 

' Three inverse boundary value problems for the heat equations in one space di- 

■ mension are considered. Those three problems are: extracting an unknown interface 

in a heat conductive material, an unknown boundary in a layered material or a ma- 
terial with a smooth heat conductivity by employing a single set of the temperature 
and heat flux on a known boundary as the observation data. Some extraction for- 
<*' , mulae of those discontinuities which suggest a relationship between the travel time 

^ I of a virtual signal and the observation data are given by applying the enclosure 

Tlj" ' method to the problems. 
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^ I space dimension. First we recall the problem. Let a > and T > 0. Let u = u{x,t) be a 
solution of the problem: 

ut = u^x in]0, a[x]0, T[, 

u^{a,t) = OfoTte]0,T[, (1.1) 
u{x,0) = in]0, a[. 

Then the problem considered therein is: extract a from a single set of the data u{0, t) and 
u^{0,t) for < t < T. 

This is a simplest one space dimension version of the problem of domain determination 
which is a typical inverse boundary value problem for the heat equation and related to 
the thermal imaging of unknown discontinuity such as cavity, defect or inclusion inside 
a heat conductive body. There are extensive studies for the uniqueness and stability 
issues of this type of problems in multi dimensions. See [1, 2] and references therein for 
several results on the issues. However, in our opinion, seeking an analytical formula that 



1 Introduction 

In [11] we considered an inverse boundary value problem for the heat equation in one 
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directly connects information about discontinuity with the data also yields another view 
for understanding of the problems. 

Recently, some new analytical methods for the such type inverse problems that the 
governing equations are elliptic were introduced. In particular, the probe method ([5, 10]) 
and factorization method {[lA, 15]) gave ways of extracting unknown discontinuity from the 
data that are given by the Dirichlet-to-Neumann map in inverse boundary value problems; 
the far field operator or the restriction of the scattered fields onto a sphere surrounding 
unknown discontinuity that are exerted by infinitely many point sources located on the 
sphere in inverse obstacle scattering problems. 

These methods require infinitely many data. By the way, the enclosure method also gives a 
way of extracting the convex hull of unknown discontinuity from the Dirichlet-to-Neuman 
map ([7]). However in some cases, using the idea of the enclosure method, one can give 
a way of extracting unknown discontinuity by a single set of the Dirichlet and Neumann 
data (see [6, 8]). In particular, the result in [8] gave a constructive proof of a uniqueness 
theorem established in [4] and the numerical implementation of a reconstruction algorithm 
of the convex hull of unknown polygonal inclusion or cavity has been done in [12, 13]. Thus 
it is quite interesting whether the enclosure method still works for the inverse problem 
for the heat equation mentioned above. In [11] we found a simple extraction formula of 
a by using the idea of the enclosure method. Needless to say, this result gives a new, 
constructive and simple proof of the uniqueness theorem: the data m(0, t) and Ux{0, t) for 
< t < T uniquely determine a (under a suitable condition on ^^.(O,^)). However, it 
should be pointed out that the result in [11] is against the past experiences in the study 
of inverse boundary value problems for elliptic equations. Here we review the result and 
point out the difference. 
Let c be an arbitrary positive number. Set 



-cr 1^1 +^^1 - -i- j, r>c-^. (1.2) 
Let 

V{x,t) (1.3) 

Given s e R we introduce the indicator function Ic{t; s): 

I At- s) = e^' r (-vJO, t)u(0, t) + ujO, t)v(0, t)) dt, r > 
Jo 

where u satisfies (1.1) and v is the function given by (1.3). 

Assume that we know a positive number M such that M > a. Let c be an arbitrary 
positive number satisfying 2Mc < T. Assume that Ux{0,t) = l(note that this is just for 
simplicity of description). Then we have the formula 

lim MMlM = .2oa. (1.4) 

r — >oo j- 

and the following statements are true: 
if s < 2ca, then lim |/c(t: s)\ = 0: 

T >00 

if s > 2ca, then lim |/c(t; s)\ = oo. 
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Why is this interesting? The reason consists of two points. 

The first point is an unexpected asymptotic behaviour of the indicator function. More 
precisely, integration by parts gives 

IJt;s) - -e^' u{a,t)v^{a,t)dt + 0{e-^^'^-'^). 
Jo 

The function e'^^v has the special character 

• if s < cx + t, then hm^ >oo e'^^\v{x, t)\ = 

• if s > cx + t, then lim,- >oo e'^'^\v{x, t)\ = oo. 

Therefore if T > ca > s, then \Ic{t;s)\ is exponentially decaying as r — > oo. Since 
e'^^\vx{a,t)\ is exponentially growing in the case when T > s > ca and < t < s — ca, 
usually we expect |/c(t; s) \ is exponentially growing as r — > oo. Then a past experience 
suggests that the right hand side of (1.4) should give —ca which has the meaning: 



-ca = sup{ . • ^ I e]a, oo[x]0, T[}. 




This is nothing but the value of the support function for the unknown domain ]a, oo[ x ]0, T[ 
at the direction (— c, —1)'^. However, in fact, we obtained —2ca which has the meaning: 



2ca = sup{ • ^ \{x,t) e]2a, oo[x]0,T[}. 




Since the u of (1.1) can be extended as a solution of the heat equation onto ]a, 2a[ by 
the reflection x i — > 2a — x, one may think that this is because of the simple Neumann 
boundary condition at x — a. However, in [11] we saw the same phenomenon to the 
case of the Robin boundary condition at x — a. So we guess that this is a universal 
phenomenon. Based on this belief, therein we gave another interpretation of 2ca. It is the 
travel time of a virtual signal with an arbitrary flxed propagation speed 1/c that starts 
at the known boundary x = and the initial time t = 0, reflects at another unknown 
boundary x — a and returns to the original boundary x — 0. 

The second point is: needless to say, the formula (1.4) yields a new, simple, construc- 
tive proof of the uniqueness theorem. 

• Let T be a fixed arbitrary positive number. For j = 1,2 let Uj satisfy (1.1) with 
a = aj{> 0) and {uj)x{0,t) = 1 for aU t e]0, T[. If ui{0,t) = U2{0,t) for aU t ejo, T[, 
then ai = a2. 

A standard and traditional proof of this type uniqueness theorem (see, e.g., [1]) is done 
by using a contradiction argument and starts with assuming, say Oi < a2- Then the 
uniqueness of the lateral Cauchy problem for the heat equation gives Ui{x, t) = U2{x, t) for 
{x,t) e]0, ai[x]0, T[. This yields that U2 satisfies {u2)x{o,i,t) = {ui)x{ai,t) = 0. Since 
U2 satisfies the heat equation in ]ai, a2[ x ]0, T[, the Neumann boundary conditions at 
X — ai, a2 and the initial condition U2{x, 0) = 0, we obtain U2{x, t) — in ]ai, a2[ x ]0, T[. 
Then the unique continuation theorem for the solution of the heat equation gives U2{x, t) — 
in ]0, a2[ x ]0, T[ and this thus yields {u2)x{0,t) = 0. Contradiction. 
Clearly this argument can relax the condition on ^^.(O,^) and tells us the importance of 
the uniqueness of the lateral Cauchy problem or the unique continuation theorem for 
the heat equation. However, this type proof gives no information about how to extract 
unknown a from the data u{0, t) and Ux{0, t). 
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It should be pointed out that, in [1] another argument in the case when T — oo and one 
has the data u{0,t) and Ux{0,t) for all t e]0, oo[, is introduced. Starting with assuming 

ai < a2, we see that U2 satisfies the Neumann boundary conditions at a; = ai, a2- This is 
same as above. However, here they do not make use of the initial condition U2{x, 0) — 0. 
Instead they have the identity 

d /■'*2 ra2 

— U2(x,t)dt^ {U2)xxix,t)dt ^0. 
Ctt J a\ Jai 

One the other hand, using the eigenfunction expansion, they show that 

lim I / U2{x,t)dt\ = 00. 

Contradiction. The proof is quite interesting, however, clearly in this argument the as- 
sumption T = cx) is essential and the same comment works also for this proof. 
Summing up, one should find another type proof that tells us the information about how 
to extract unknown a from the data u{0,t) and u^,{Ojt) for t g]0, T[ with T < cxo. We 
think that our proof presented in this paper gives an answer to this natural question. 

The aim of this paper is: to confirm further that the interpretation in the first point 
still works, at least, in one space dimensional case by considering three typical inverse 
boundary value problems for the heat equations. Those three problems are: extracting an 
unknown interface in a conductive material, an unknown boundary in a layered material 
or a material with a smooth conductivity. 

In a future study we will consider the multidimensional version of those problems. 
Remcirk 1.1. In this paper we always consider the solutions of the heat equations in the 
context of a variational formulation. In particular, every solutions in this paper belongs 
to the space W{0, T; H^ fl), (H^n))') with m^(0, t), M^(a, t) G L\0,T) where Vl =]0, a[ 
and satisfies the governing equation in a weak sense. We refer the reader to [3] for the 
detail. 

2 Statement of the results 

2.1. Extracting interface 

Let < b < a. Define 

{71, if < X < b, 
72, lib < X < a 

where both 71 and 72 are positive constants and satisfies 72 7^ 71. 
Let u be an arbitrary solution of the problem: 

ut = (lUx)x in]0, a[x]0, T[, 

(2.1) 

u{x,0) = Oin ]0, a[. 

We assume that both 71 is known and that a, b and 72 are all unknown. 
Inverse Problem A. Extract b from u{0, t) and 71 Ux{0, t) ior < t < T. 
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Let c be an arbitrary positive number. Let 

v{x,t)^e-^^^<i/{x,z) (2.2) 
where ^{x, z) — e^^% Zi — zj ^fy^ and z is given by (1.2). 

The function t> is a complex valued function and satisfies the backward heat equation 

"Vt + 7i Vj^j. = in the whole space-time. 

Definition 2.1. Given c > define the indicator junction Icij) by the formula 



)= r(-7i^^a:(0,t)M(0,t) + 7iM^(0,t)^;(0,t))(it, r>c- 
■/o 



where u satisfies (2.1) and v is the function given by (2.2). 
Define 



w 

Then this w satisfies 



{x) — w{x^T^— / e u(x,t)dt, < X < a. 
Jo 

{-fw'y -z''w^e-''''u{x,T) in]0, a[ (2.3) 



and 

w 



'(0)= / e-'^%^{0,t)dt, w'{a)^ [ e-'^%^{a,t)dt. 
Jo Jo 

Our first result is the following theorem: 

Theorem 2.1. Assume that we know a positive number M such that M > 2b/ ^y^. Let 
c be an arbitrary positive number satisfying Mc < T. Assume that 

,w'(a), ( , b {a — b)\^ ,^ 
lim — H- exp cr ^ ^ =0 2.4 

and there exist a positive constant C , positive number Tq, real number /i such that, for all 

T > To 

Ct" < \w'{0)\. (2.5) 



Then the formula 

lim i^S^ = -2 ^, 

is valid. 

Note that the boundary condition at a; = a is not specified. However, the condition (2.4) 
implicitly restricts a possible boundary condition at a; = a. The situation dramatically 
changes in the case when (a) b/ ^/ti < {a — b)/ ^/72 or (b) b/ ^/ti > {a — b)/ ^J^. 
In case (a) the condition (2.5) automatically ensures that (2.4) is vahd since we have always 
u''(a) = 0(1) as r — > oo(Remark 1.1). This is reasonable under our interpretation since 
the signal started from a; = a at t = arrives at a; = after the arrival of the signal 
started at a; = at t = (see Figure 1). 

However in case (b) to ensure (2.4) w'{a) has to decay exponentially as r — > oo. This is 
a strong restriction and can be interpreted as a condition that kills a signal started from 
a; = a at t = (see Figure 2). 
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Figure 1: (a) < (a - h)/^. 




The condition (2.5) gives a restriction on the behaviour of the flux Ux{0,t) as t | 0. 
Let S satisfy < S <T and m > be an integer. Since 



^^ = T + ^2cV^/l-^, 

V C^T 



a change of variable yields, as r — > oo 



1 i-rs -^(l + i2cVWl- ^ 



/ e-' H"'dt = —— e V cV 



1 ^.oo -cMi + n/i-^)'^ 



oo 



^2(m+l) 

Integration by parts gives 



1 



/ e V V^dt^mlKiT 

Jo 



m+l 



where 

Therefore we obtain, as r — oo 

rS 



^2(m+i) /■ e-^^t^m^^ _^ rn\{-^r+^ + 0. 

t/ 



It is easy to see that if / e -^>^(0, T) and for a suitable positive constant C, < Cf^' 

a. e. in ]0, 5[, then, as r — > oo 



JO 

Now assume that, we have, for some 5, a positive constant C 

|ii,(0, i) - (^^i' + + • • • + < Cr'+\ a. e. in ]0, 5[ 

where n, / and m! are integers and and satisfy n > / > and m' > 2/ + 1; ... g(„ are 
constants and satisfy gi ^ 0. Then from the computation above we see that, as r — > oo 



^2(i+l) 
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This implies that Ux{0,t) satisfies (2.5). In particular, the condition (2.5) is satisfied if 
Ux{0,t) satisfies one of the following conditions 

• u^{0,t) e C2([0, 5]) and u^{0,0) ^ 

• Ux{0, t) G C^'+^([0, 6]) with / > 1 and t = is the zero point of Ux{0, t) with order /. 
2.2. Extracting unknown boundary. Layered material 

Let = 6o < &i < &2 < • • • < 



7(x) 



a. 



7i, 



72, 



iibo < X < bi, 
ifbi < X < 62, 



. 7m, if bm-l <X <br, 



where 71, 72, • • • , 7m are positive constants and m > 1. 

Let a > 0. Let u — u{x, t) be an arbitrary solution of the problem: 

ut = i'JUx)x in]0, a[x]0, T[, 

-fmUx{a,t) + pu{a,t) = fort e]0, T[, 

u{x, 0) = in]0, a[ 



(2.6) 



here p > is an arbitrary fixed constant. 

Inverse Problem B. Assume that 71, ■ ■ -, 7„j are known and that both of p and a are 
unknown. Extract a from ^(0, t) and 71 Ux{Oj t) for < t < T. 

In this subsection, instead of in the previous subsection the function ^' in the 
following which is a consequence of Lemma 5.1 in Section 5 plays the same role; z is given 
by (1.2) and set zj — z/^/^, j — 1, • ■ ■ ,111. 

Proposition 2.2. There exists a positive number C — C{c, 71, • • • , 7m) (> c~^) such that: 
for each t > C there exist a unique U — {Bi, A2, B2, ■ ■ • , ^m-i, -^m-i, ^m) such that the 
function ^ defined by the formula 



{ e^^i + Bie- 



*(x) = *(x;c,7i,---,7m,r) = < 



ifx< bi, 
ifbi<x< 62; 



satisfies the equation 



Am-ie'''"'-' + Bm-ie ""^^-^ if bm-2 <x < bm-l, 

, ^m^^^"*, 'if bm-l < X 

(7*')' - -z^* = mR 
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where 



j{x) = < 



71, if-oo<x<bi, 

72, ifbi<x< 62, 



Using "if in Proposition 2.2, we define tiie special solution of the backward heat equa- 
tion vt + {'yvx)x = in R by the formula 

v{x, t) = e-^'**(x; c, 7i, • • • , 7^, r), r > C{c; 71, • • • , 7,„) 

Now we can define the indicator function. 

Definition 2.2. Given c > define the indicator function Ic{t) by the formula 

W = r (-^1 ^^(0' ^) + ^1 ^^(0' ^)) ^ > • • • ' ^"^) 

where u satisfies (2.6). 

The following gives a solution to the problem mentioned above and generalizes Theo- 
rem 4.1 of [11]. 

Theorem 2.3. Assume that we know a positive number M such that 



M>2 



bi 



m— 1 



(6, 



'i+i - "i. 



71 ^ § 



Let c be an arbitrary positive number satisfying Mc < T . Assume that there exist positive 
a constant C, positive number tq{> c"^), real number /i such that, for all t > tq 



Then the formula 



lim 

T >0O 



Cr'' < I r e-'''uJO,t)dt\. 
Jo 



loglWI _ o J fei I y jbj+i-bj ) 



(2.7) 



(2.8) 



is valid. 

Needless to say, the quantity 

.J^^j^ib^) 

can be interpreted as the travel time of a virtual signal with propagation speeds 
\/7i/c, \/72/c, ■ ■ ■ , y/Tm/c in the layeres < a; < 61, 61 < x < b2, ■ ■ ■ , &m-i < x < a, 
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respectively that starts at the known boundary x — and the initial time t — 0, reflects 
at another unknown boundary x — a and returns to the original boundary a; = 0. 
2.3. Extracting unknown boundary. Material with smooth conductivity 

Let a > 0. Let M > a. Let 7 E ^^([O, M]) and satisfy 7(x) > for all x e [0, M]. 
Let u = u{x, t) be an arbitrary solution of the problem: 

ut = {-fu^)x in]0, a[x]0, r[, 

7(a) u^{a, t) + pu{a, t) = for t e ]0, T[, (2.9) 

u{x, 0) = in]0, a[ 
here p > is an arbitrary fixed constant. 

Inverse Problem C. Assume that both M and 7 are known and that both of p and a 
are unknown. Extract a from u{0, t) and 7(0) Ux{0, t) ior < t < T. 

We start with a fact which can be deduced from a combination of Theorem 1 of p. 48 
in [16] and the Liouville transform. 

Proposition 2.4. Given z with Rez < there exists a solution = ^( • ;z,M) of the 
equation {jy')' — z'^y — 0, < x < M such that, as \z\ — > 00 



<i/{x;z,M) = {7(x)}-^/^ exp ^z 



\z\ 



■^'{x-z,M)=z{^{x)}-^'Uxp \^ £^l=j + 

uniformly in x G [0, M]. 

Using \E' in Proposition 2.4, we define the special solution of the backward heat equa- 
tion vt + {'jVx)x = in ]0, M[x ]0, T[ by the formula 

v{x,t) = e"^'**(a;;2;, M). 

Now we can define the indicator function. 

Definition 2.3. Define the indicator function I{z) by the formula 

I(z)^ r (--f(0)vJO,t)u(0,t)+-f(0)uJO,t)v(0,t))dt, Rez<0 
Jo 

where u satisfies (2.9). 

The following gives two solutions to the problem mentioned above and generalizes 
Theorem 4.1 of [11]. 

Theorem 2.5. Assume that we know a positive number M such that M > a. 
(1) Let c he an arbitrary positive number satisfying 



i-M fjq- 

T>2c (2.10) 
^0 
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Let z be the number given by (1-2). Assume that there exist positive constant C and a 
positive number Tq{> C^), real number /i such that, for all r > Tq 

Ct" < I r e-'"'uJO,t)dt\. (2.11) 

Then the formula 

lim ^osmi ^ (2.12) 
r Jo ^ 

is valid. 

(2) Assume that there exist positive constant C and a positive number tq, real number 
II such that, for all r > tq 

Cr^ < I r e-^''uJO,t)dt\. (2.13) 

^0 



Then the formula 

is valid. 
The quantity 



lim i^g^(zl^^_2 r_^^, (2.14) 



„ dx 
ic 



■f 

Jo 



coincides with the travel time of a virtual signal with variable propagation speed \J^[x)/c 
that starts at the known boundary x = and the initial time t = 0, reflects at another 
unknown boundary x = a and returns to the original boundary a; = 0. 

The condition (2.13) is less restrictive compared with condition (2.11). We can easily 
see that the condition (2.13) is satisfied if ^^.(O, t) satisfies one of the following conditions 

• for a positive constant C Ux{0,t) > C a.e. in ]0, S[ 

• u^iO,t) e CH[0, 6]) and uj{0,0) ^ 

• Ux{0, t) e C'"'"^([0, 6]) with / > 1 and t = is the zero point of Ux{0, t) with order /. 
See the end of subsection 2.1 for the comparison. 

As a corollary of Theorem 2.5 we obtain a direct proof of uniqueness theorem: the 
data u{0,t) and Ux{0,t) for < t < T uniquenely determine a provided one of (2.11) for 
a c satisfying (2.10) or (2.13) is satisfied and both M and 7 are known. Note that T is an 
arbitrary fixed positive number. In the proof we never make use of the uniqueness of the 
lateral Cauchy problem nor the unique continuation theorem for the heat equation with 
a variable coefficient. 

3 Proof of Theorem 2.1 

Since ^ satisfies 71 — 2;^^ = in R, integration by parts gives the expression of the 
indicator function: 

Ic{^) = -71 ^'{a)w{a) + 72 w'(a)*(a) + (71 - 72) / w'{x)-^'{x)dx 

Jb 

(3.1) 

-e-^^ ru(e,r)^(e)de. 

Jo 
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Let y be the solution of the boundary value problem 

(7y')'-^'z/ = 0in]0, a[, 

y'iO)^w'iO), y'ia)^w'ia). 
Define the principle part of the indicator function 

7°(t) = -7i*'(a)y(a) + 72Z/'(a)*(a) + (71 - 72) T y'ix)^'{x)dx. (3.2) 

It is easy to sec that Theorem 2.1 is a direct consequence of the following two lemmas 

Lemma 3.1. Assume that w'{a) and w'{0) satisfy (2.4) and (2.5). Then, choosing 
suitable positive constants C[ and C'2 and positive number Tq, we have 

C'^t" < |/°(T)|e2'=^''/V7r < c'^^ Vr > r^. 
Lemma 3.2. As r — > 00 

W = /°(r)+0(e-^^). 

First we prove Lemma 3.1. The y has the expression 

Ae^^i + Sie-^^\ iiO<x<b, 



Aoe'^'^ + Boe-^'^ ii b < x < a 



y{x) = y{x,T) 

'i2 e - -D2 ' 

and has to satisfy the transmission conditions: 

y{b - 0) = y{b + 0), 71 y'{h - 0) = 72 y\h + 0). 
These and the relation Z2/Z1 — ^71/72 yield the system of equations: 

Aie^^i +5ie-^"i = ^2 6^^^ + 52 6-^"^ (3.3) 

Ai e^^i - Si e-^'^i = (^2 e''^^ - S2 e-''^^). (3.4) 

V 7i 

Moreover the boundary conditions y'(0) = w'{Q) and y'(a) = w'{a) yield 

^1 (Ai - Si) = «;'(0) (3.5) 

and 

Z2 {A-2 e"'^ - B2 = w'{a). (3.6) 

A combination of (3.3) and (3.4) gives 
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(3.7) 



Define 



We have 



'-kl 



, R, 



kl 



Ik - 

ik + ^/Yi' 



Then (3.7) becomes equations 



1 

T\2 



-t 12 



-^12 4 hz2 , 1 

77^ ^2 6 + — 
-'12 -tl2 

Substituting (3.8) and (3.9) into (3.5) and (3.6), we obtain 



R 







V B2 e-''^ ) 





I2W {a} 



where 



R 



( g(a-6)22 _g-(a-&)z2 ^ 

/I Rl2\ 

(1 



— e 



26 zi 



) 



V 

/ g{«-&)22 



T12 




R\2 


1 







V ^{1-Ri2e'''^) -^(i?i2-e^^- 

\ -^12 -'12 



_g-(a-6)z2 \ 



A direct computation gives 

(deti?) ri2 e^"-'')^^ = 1 + i?i2(e2^"-'')^2 - e^*"^!) - e2''^i+2(«-'')^2 
This yields, as r — 00 

(det R) T12 e("-^)^^ = 1 + (e-2c^ (VVtT, (a - 6)/V^) 



(3.8) 
(3.9) 



(3.10) 



(3.11) 



Therefore R is invertible for sufficiently large r. It follows from (3.10) that A2, B2 have 
the form 



-bz2 



Bo = 



^(deti2) Vri2 



— (i?i2 - e^^^i) v^w'(a) + e-(«-^)^^ V7rw'(0) e^^^ 



^(deti?) Vri2 



— iR^2 e''^^ - 1) ^2w'{a) + e^"-^)^^ V7r^/;'(0) e'^' 



(3.12) 
(3.13) 
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Using (3.12) and (3.13), we get two crucial formulae: 

(71 - 72) [y'{xW{x)dx = {C V^^'(a) + DT,^ Vll^'m 

where 

C = (VtT - ^){Ri2 - e'"'') (e«^i+(«-^)^^ - e'^') 
+{VYi + VY2){Ri2e'''' - 1) (e»^i-(»-^)^^ - e^^i), 
D = 2^e("+^)^^ 

-7i^''(a)2/(a) + ^2y'{a)^{a) - ^^^^^^^^^ (C'V7^w'(a) + D T,2 V^w' 

where 

C = -{VYi-VT2){Ri2 - e2^^i)e'^^i+(«-^)^^ 
-(V7r+A/^)(^i2e'^^i - l)e-l-(«-^')^^ 

A combination of those and (3.2) gives the representation formula 
(deti?)Ti2/°(r) = -2^,^w'(a) e'^' 

- ((v^ - v^) e^''^^-^"-")^^ + (^71+ v^) e2''-i+(«-^)^^) T^^ V^iw\0). 
Note that 



(3.14) 



-2cr (a— 6)/y^ 



Note also that Remark 1.1 gives w'(0) = 0(1) as r — > 00. Using these, (3.11), (3.14), 

(2.4), (2.5) and the assumption 71 7^ 72, we get the assertion of Lemma 3.1. 

□ 

Next we give a proof of Lemma 3.2. Recalling (3.1) and (3.2), we get 

/c(r)-/°(r) 

= -71 *'(a)e(a) + (71 - 72) / e'ix) ^'ix)dx - e"^ ^ / ui^, T) ^(0^^ 

Jb Jo 

where e{y) — w{y) — y{x). One knows 

/\(e,r)*(0 = o(i). (3.16) 

The e satisfies 

(7 e')' -z^e^ e-^' ^^(a;, T) in ]0, a[, 
e'(0) = e'(a) = 0. 
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Multiplying the both sides of the equation by e and integrating the resultant over interval 
]0, a[, we have 

Jo Jo Jo 

Since 



9 ' 



a standard argument yields, as r — > oo 

||e|ki(]o,a[) = 0{e-^^). 
A combination of this and the embedding i?^(]0, a[) C C°([0, a\) gives the estimates 

-7i*'(a)e(a) = 0(Texp (-r(r + ^)V, 



Hi 

(3.17) 

r e'(a;) ^'(x) dx = 0(t exp ( -r (T + 4^) ) ). 



. 7i . 

Now we obtain from (3.15), (3.16) and (3.17) the assertion of Lemma 3.2. 
□ 



4 Proof of Theorem 2.3. Part 1. Asymptotic be- 
haviour of w{a) 

Define 

wix) — w(x, t) — I u(x, t) e ^ * dt, < X < a. 
Jo 

This w{x) satisfies 

(7^7 -z'^w^ e-^'^ii(x,T) in]0, a[, 
^rnw'{a) + pw{a) = 0. 

In this section we study the asymptotic behaviour of w{a) as r — > oo. For the purpose it 
suffices to study the asymptotic behaviour of the solution of the boundary value problem 
as r — > oo: 

(7y')'-^'z/ = 0in]0, a[, 
y'{0) = w'{0), 

imy\a) + py{a) = 0. 

This is because of 

Lemma 4.1. Let p>0. The formula 

w(a)^y(a) + 0(e-^^), 
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is valid. 
Proof. Define 

e{x) — w{x) — y{x), < x < a. 

This function satisfies 

(7 e')' -z'^e = e"^" ^u{x, T) in ]0, a[, 
e'(0) = 0, 
7^e'(a) +pe(a) = 0. 

Hereafter a combination of a standard argument and the embedding if ^(]0, a[) C C°([0, a]) 
yields the desired estimate. 

□ 

For each j — 1, • • • , m one can write 

y{x) = Ajc'''^' + Bje~''^', bj^i < x < bj. 
From the equation, it follows that 

y{bj-0)^yj+i(bj + 0), 

7il/;(&i-0)=7,+il/;+i(&i + 0). 
Therefore the coefficients Ai, Si, • • • , Am, have to satisfy the system of equations: 

zMi-Bi)=w'{0); (4.1) 

for each j = 1, • • • , m — 1 

AjC^^'^ + Bje'^^^^ = Aj+ie^^^^+^ + Bj+ie'^^'^+\ 

^Mj^''^' - Bje-"^'^) = 7,-+i ZJ+l{A^+^e''^^^+^ - Bj+.e-'^^^^^); 
7^ ZmiAme'''"^ - S^e-"^-) + plA^e"^- + S^e-"^-) = 0. (4.3) 



(4.2) 



Set 







1 








\ 


Bje-^^'^ 1 


V 



Hi 



1 \ 

e-^(''^-i-^^) Zj ) 



Then one can rewrite the equations (4.1), (4.2) and (4.3) in the matrix form: 

(1 -e^^- ) Xi = ^^^^^^e^-^ (4.4) 
^ ' z 

e('^-'^+^)^^+^ir,+ia,+iX,+i = K,X,, (4.5) 

((V7m + ^) -(V7m-^) )x„ = 0. (4.6) 
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Set 

L{z) = (Xf' K, a2) • • • {K-\ aj. 

Prom (4.5) we have 

m—l 

J2 ipj - bj+i) Zj+i 
Substituting this into (4.4), we obtain 



m—l 



Zl„ j = l 



(4.7) 



The equations (4.6) and (4.7) are equivalent to the equation 



m—l 



Xm — 



y ( 1 -e^^'i^i ) L{z) J 



(4.8) 



Solving equation (4.8), we obtain Xm. The problem is the asymptotic behaviour of L{z) 
as r — > oo. 

Define the transmission coefficient T^i and reflection coefficient R^i by the formula 



- fcZ 



/7ik + 

Set 5 = mmj=i^....rn{{bj - bj-i)/^}{> 0). 
Lemma 4.2. VFe /iave, as r — >• cxo 



, R, 



kl 



Tl2 ■ ■ ■ Tjn-1, 



( 1 



V ^12 



+ 0(e-2'=^^) 



(4.9) 



Proof. Using the expression 



1 







and 



we have 



1 (9(g-2cT(6j-63_i)/^ 
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This gives 



Since 



I 1 



1 0^ 



we obtain 



yi - V72 



\ 







This is nothing but (4.9). 
□ 

As a direct consequence of (4.9) we have 



/ 



V w 

1 ^ 



/7m + 



2; 



1,(1 -e^^-^)L(z); 

/ • • • Tjn-l,m \ 



— 1 Ti2 ■ ■ ■ Tm-1, 

and (4.8) therefore yields 



V 



+ 



2p 



( ° 

V —1 Ti2 ■ ■ ■ Tjn-1,; 







o 



m—1 



Xr, 



Tl2 ■ ■ ■ Tm-l,m{ 



2p 



+ 



as r — > oo. 
Define 
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and 



X Zi, 

XZ2 + ip2, 



if — oo < X < bi, 
iibi < X < 62, 



, XZm + ^m, if bm-1 < X < OO. 

The function (p has the unique continuous extension to the whole real line since (f{bj — 0) 
ipipj + 0) for each j = l,---,m — 1. We denote the extension by </? again. 
Since y{a) = + (X^)2, a = &m and 

m—l 

biZi-Y, {bj - bj+i) Zj+i 



= 61 2^1 - (61 - 62) ^2 - (&2 - ^3) ^3 (bm-1 - bm) 

= bi{zi - Z2) + 62 (2:2 - ^3) H h bm-1 {Zm-1 - Zm) + bm Zm = ^{bm), 

we obtain the formula 

y^a) = e'^^'^) T12 • • • Tm-l,m{^ 

z 

Note that 

(^(a) = —z 



ZJl-n 



b, , - 6,. 



7i ^ 



(4.10) 



(4.11) 



5 Proof of Theorem 2.3. Part 2. Asymptotic be- 
haviour of ^(a) 

Integration by parts yields that the function \E' is the (weak) solution of the equation 
(7 ?/')' ~ ^'^y = in R if and only if ^ satisfies, for each j = 1, • • • , m — 1 

*(6,-0) = *(6, + 0), 

7,- *'(6, -0)=7,+i *'(&,• + 0). 
This yields the system of equations for Si, >l2, -B2, • • • , Bm-i, Am- 

g(5,-5,+,)^.+i;^.^^«.^^y.^^ = KjYj, j = 1, ■ • •,m - 1 (5.1) 

where 













/ 1 \ 


Yi^ 






V — A pbmZm 
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Prom (5.2) one has 



m— 1 



(5.3) 



J2 {bj - bj+i) Zj+i 

m—l 

^ Ame^"'^"'e^=^ L{z) 

Vo 

Set 5 = mmj=i^...^rn{{bj - bj-i)/^}{> 0). 

The following indicates the meaning of the coefficients T^i and Rki- 

Lemma 5.1. For sufficiently large r the system of equations (5.1) is uniquely solvable 
and, as r — > oo the formulae 

Aj e-^^ ^1 e«=i = • • • r^--i, ^ + 0{e-'^''^^) (5.4) 

anc? 

i-i 

Y^ipi - bi+i) zi+i 

Bje~''''e-'^'^ei=i ^ Rj^+iT^- ■ -Tj.^j + 0{e-^''^^), (5.5) 

are valid. 

Proof. Using (4.9) and (5.3), we find for a sufficient large r and obtain (5.4) for j = m. 
Next (5.1) for J = m — 1 yields that Am-i and Bm-i are uniquely determined and that 
the formulae (5.4) and (5.5) for j = m — 1. For general j we make use of the recurrence 
formulae 

which are equivalent to (5.1). Note that we set A^ — l and Bj^ — 0. 

□ 
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Prom Lemma 5.1 one has, as r — > oo 

^2 6-^2 = Ti2 + 0(e-2^^"), 

Ase-'^'^T,2T23 + 0{e-'"^^), 

e-^- = ri2 r23 • • • r„_i,„ + o(e-2<=^-); 

S2 e-'^^^^ e-'^^ = ri2 i?23 + 0(e-'^^^), 



e-'^-i = ri2 T23 • ■ ■Tm-2,^-lRm-l,m + 0(6-2^''^). 

Moreover, it follows that 



This gives, as r — oo the asymptotic formula of ^: 



Ti2e^(^), 



if — oo < a; < 61, 
if 61 < X < 62, 



, 7i2 • • • r„_i,,„ 6"^^^^ if < X < 00. 
The formula (5.4) for j = m gives 

vl>(a) e-'^('') = Ti2 ■ ■ ■T„_i,„ + ©(e^^'^''^). 
The following estimates are a direct corollary of (5.4) and (5.5): 

/ 



1^,1 = 



j-1 ^^ 



J 



\Bj\ = 



V 



(5.6) 
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Applying these estimates to the expression 

Jbj—i Jbj—i Jbj—\ 

we obtain the estimate ^ 

/ u{i,T)^{m = 0{l). (5.7) 

RecaUing the definition of ty, the equation ^'a) — Zm^{a) and using integration by 
parts we have 

Ic{r) = 71 «;'(0)*(0) - w(0) 7i*'(0) 

Jo 

Then Lemma 4.1, (4.10), (5.6) and (5.7) yield the asymptotic formula for the indicator 
function. 

Proposition 5.2. As r — oo the formula 

Ic{r) e-'^^^^ = -2 VlPuyj'iO) (ri2 • • • Trn-i,mni + + 0(1)} 

-r{r-2c(-^+ V^^±i— ^)} 
+0(re ^ ), 

Theorem 2.3 is an immediate corollary of (4.11), Proposition 5.2, the expression 



w'(0) = / e-''^%JO,t)dt, 
Jo 



the assumption (2.7) and the fact ^^'(O) = 0(1) as r — > oo. 

Remark 5.1. Once (p{a) is known (see (4.11)), then one immediately obtains the extrac- 
tion formula of p: 

( / (r) '^(''^ \ 

6 Proof of Theorem 2.5 

Define 



w 

This w satisfies 



(x,t) = / *u(x,t)dt, < X < a. 
Jo 

{-fwy-z'^w^e'-^^^u{x,T), in]0, a[, 
7(a) i(;'(a) + pw{a) — 0. 
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Integration by parts gives the expression 

I{z) = -w{a,z){-f{a)^'{a;z,M)+p^{a;z,M)) - T e-^^^u{x,T)^/{x; z, M)dx. (6.1) 
Proposition 2.4 gives 

7(a) ^\a; z, M) + p^{a- z, M) = -^e^^'" ^^{1 + O j }. (6.2) 

Therefore it suffices to study the asymptotic behaviour oiw{a,z) as \z\ — oo. For the 
purpose first we study the asymptotic behaviour of the unique solution of the boundary 
value problem as — > oo: 

(7l/0'-^'?/ = Oin]O, a[, 

y\o) = 1, 

7(a)|/'(a) + py{a) = 0. 
We make use of the Liouville transform. Define 

1 dx 



and 

Ka 

Denote by a;(s; a), < s < it the inverse of the function s = s{x; a). Set 

y{s) = X„7(0)^/%(x(.; a)y/'y{x{s; a)). (6.3) 



, , 1 dx 
six: a) = —— \ I , \) < X < a. 



Then this y satisfies 



where 



y"-(K„V + ^„(s))y = Oin]0, 7r[, 
y'(0)-/i„y(0) = l, 
y'(7r)+^ay(7r) = 

7'(0) 4p-7'(a) 



and 



4i^„ 7(0)3/2' 4i^„7(a)3/2' 

/„(.)=7(^(«;a))'/', ^a(«) = f^- 
Lemma 6.1. ^5 \z\ — >• oo we have 



KaZ 



y(7r) = ^e^--{l + o(^)}. (6.4) 
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Proof. It is easy to see that y has the expression 

vis) 



where (f is the unique solution of the initial value problem: 

ip" -{Klz^ + g{s))^ = Om]0, 7r[, 

ip{7r) = 1, 
<^'(7r) = -H,. 

Here we cite a known important fact: there exists a fundamental system of solutions 
ei(x, z) and 62(0;, z) of equation y" — {K^ z^ + ga{s))y = in ]0, tt[ such that as \z\ — > 00, 
Kez < 0, uniformly in x, 

e,{s,z)=e'''^''{l + o(^)}, 



e[{s,z 



e2(.,^) = e-^''-{l + 0(^^j}, 
e',{s,z) = -Kaze-'^'^'^I + O^-^^} 



See again Theorem 1 of p. 48 in [16]. 
Thus one can write 

ip{s) ^ Ai ei{s, z) + A2 62(5, z) 
where Ai and A2 are constants. The initial conditions on 93 yields 

iy(ei(-,z),e2(-,^))(7r) 

W{ei{- ,z),e2{- ,z)){t:) 
The asymptotic behaviour of ei and 62 yields 



W{e^{ ■,z),e2(-, z)){7r) = -2X„ z{l + O ( ^] } 



\z\ 



and also 
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Prom these we obtain 



\z\ 



z\ 



This yields (6.4). 

□ 

A combination of (6.3) and (6.4) gives 



Using a standard argument, we have 



^^"'"^ ■ 0(e-R^^'^). (6.6) 



w'{0, z) 

A combination of (6.5) and (6.6) yields 

•"(". ^) = Kl'^lultay,: ^^^ - O (r)} + 0(e-^^^'T^X,, .)). (6.7) 
Now from (6.1), (6.2) and (6.7) we obtain the crucial formula: 



\z\ 



+0 (e-^^ ^|^'(0, z)\\z\ + O [e'^^ ^) 



(6.8) 



Then Theorem 2.5 follows from this formula (6.8), for the choices z — —ct (1+i^l - l/cV) 
(case (a)) and z — — r (case (b)) and the expression 



w'(0,2) = / *M^(0,t)cit. 
h 

More precisely, both (2.11) in case (a) and (2.13) in case (b) ensure 

Ct^ < \w'{0,z)\, T>T^. 

Prom Remark 1.1 one has w'{Q,z) = 0(1) as r — > oo. These yields the estimate: for 
suitable positive constants Ci, C2 and for all r >> 1 

Cit" < |/(^)|e-2^aRez7r < 
Note that (2.10) is essential in the case (a). This completes the proof of Theorem 2.5. 
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